Energy dependent growth of the nucleon and hydrodynamic initial conditions 
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Due to gluon saturation, the growth of the inelastic nucleon-nucleon cross section with increasing 
colhsion energy ^/s results in a broadening of the nucleon's density distribution in position space. 
This leads to a natural smoothing of the initial energy density distribution in the transverse plane of 
the matter created near midrapidity in heavy-ion collisions. We study this effect for fluctuating initial 
conditions generated with the Monte Carlo Kharzeev-Levin-Nardi (MC-KLN) model for Au+Au 
collisions at the Relativistic Heavy Ion Collider (RHIC) and the Large Hadron Collider (LHC). We 
argue that at the LHC viscous hydrodynamics is applicable at earlier times than at RHIC, not only 
because of the higher temperature but also since the length scale over which the initial pressure 
fluctuates increases with collision energy. 

PACS numbers: 25.75.-q, 12.38.Mh, 25.75.Ld, 24.10.Nz 



I. INTRODUCTION 

Hydrodynamics has been successful in modeling 
the hot, dense matter produced in ultra-relativistic 
heavy-ion collisions, particularly in describing the 
large collective flow observed at RHIC and the LHC 
[1]. A fully quantitative description of the experi- 
mental data requires, among other ingredients, real- 
istic initial conditions. These initial conditions cur- 
rently remain a significant source of uncertainty in 
predicting final state obervables. 

Optical models [2-4] for the 2-dimensional ini- 
tial density profile in the plane transverse to the 
beam axis, generated from the overlap of longitu- 
dinally integrated Woods-Saxon distributions ("nu- 
clear thickness functions" ) , yield radially symmetric 
entropy and energy densities for zero impact param- 
eter (b = 0) collisions. This symmetry drives radi- 
ally symmetric hydrodynamic expansion, with zero 
anisotropic flow coefficients 7j„ = (cos(n0)) (where 
is the azimuthal angle in the transverse plane rel- 
ative to the impact parameter, and the average is 
taken with the azimuthal distribution of the final 
state particles). The non- vanishing elliptic flow 
observed in central (6 ~ 0) Cu-|-Cu and Au-|-Au colli- 
sions at RHIC [5] and the non-zero harmonic flow co- 
efficients wi , . . . , We measured in central Pb+Pb col- 
lisions at the LHC [6-8] are in direct confiict with 
this prediction and require new physics to model the 
density of matter produced in relativistic heavy-ion 
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collisions. 

Collective flow is the hydrodynamical response 
to the pressure gradients in an individual collision 
event, and the latter fluctuate from event to event. 
In essence, each nuclear collision makes a quantum 
mechanical measurement of the actual transverse nu- 
cleon positions inside each nucleus at the time of 
collision, and the outcome of this measurement fluc- 
tuates stochastically even though the corresponding 
probability distribution (the average nuclear density 
profile) is smooth. Therefore, the initial density dis- 
tribution of matter produced even in a perfectly cen- 
tral (6 = 0) collision is not azimuthally symmetric 
(and the resulting anisotropic fiow coefficients w„ do 
not vanish) even though the ensemble-averaged ini- 
tial density has this symmetry. 

This was first pointed out by Miller and Snellings 
[9] who demonstrated the importance of event-by- 
evcnt initial state fluctuations on anisotropic collec- 
tive fiow using a Monte Carlo implementation of the 
Glauber model [10]. A Monte Carlo version of the 
KLN model [4] , which incorporates the idea of gluon 
saturation at high collision energies, was later de- 
veloped by Drescher and Nara [11] and further im- 
proved in Refs. [12-14]. 

The original MC-KLN code [11, 12] modeled nu- 
cleons as homogeneous cylinders along the beam di- 
rection whose cross section is determined by the (en- 
ergy dependent) inelastic nucleon-nucleon scatter- 
ing cross section cr^N ■ This leads to step- function- 
like discontinuities in the nuclear thickness function 
T4(r^), which (through the gluon saturation mo- 
mentum Qs(r_L) [4]) controls the initial transverse 
entropy density profile s('rj^) of the collision fire- 



ball [15]. These artificial discontinuities generate 
big fluctuations in the initial hydrodynamic veloc- 
ity gradients which create such large viscous pres- 
sure gradients that viscous hydrodynamics breaks 
down. This can be avoided by using a more realis- 
tic Gaussian profile for the nucleon density [14]. In 
[14] it was assumed that small- a:: evolution is local in 
impact parameter space, which may underestimate 
the broadening of the dipole scattering amplitude 
in impact parameter space with increasing collision 
energy. We here point out that, if one follows the 
idea of gluon saturation, the width of the gluon dis- 
tribution inside a nucleon should grow with energy, 
due to transverse diffusion of small- a; gluons: as -^/s 
increases, both and a should become larger. 

We model this by an energy-dependent broadening 
of the Gaussian nucleon thickness function Tp{r±)} 
We show that this broadening leads to a smooth- 
ing of fluctuations in the initial density distribution 
of higher-energy heavy-ion collisions and to a cor- 
responding decrease of the higher order eccentric- 
ity coefflcients e„>2 (which drive the higher order 
anisotropic flow coefficients w,i>2 [19]) at the LHC 
compared to RHIC. 



II. SIZING THE NUCLEON 

Wc consider Gaussian nucleons with the normal- 
ized density distribution (the subscript p stands for 
" proton") 



(27rS)3/2 ' 



{Jd^rp,{r)^l), (1) 



corresponding to a 3-dimensional rms radius of the 
nucleon of r™'^ ^ yHj^ = \/3B . Eq. (1) yields the 
nucleon thickness function (defined as the integral of 
the nucleon density along the beam axis) 



dzpp(r±,z) 



2ttB '■ 



I 



d^r^Tpir,^) = 1, 



(2) 



^ The possibility of measuring this broadening at a fu- 
ture electron-ion collider has recently been pointed out in 
Ref. [16]. We also note that similar ideas precede the devel- 
opment of the concept of gluon saturation and follow from 
Froissart's work on the unitarization of scattering ampli- 
tudes [17], sec [18] for a review. 



corresponding to a transverse rms radius 
r™'' = \/2B. We interpret it as the longitudi- 
nally integrated ("transverse") density of colored 
partons ( "gluons" ) inside the proton. 

In this section, we determine the Gaussian width 
i? as a function of y/s by relating it to the y/s- 
dependent inelastic nucleon-nucleon cross section 
(j'^^. IVIotivated by the Glauber- IMueller formula for 
the impact parameter dependent cross section of a 
color dipole with the quarks and gluons inside a pro- 
ton [20, 21], we make a Glauber-like ansatz [14] for 
the probability P{b) of a nucleon-nucleon collision 
at impact parameter b: 



P{b) ^ 1 - exp[-aggTpp{b)]. 



(3) 



Here Tpp{b) is the transverse density for binary col- 
lisions between gluons, 

Tpp{b) = J d\^Tp{r^)Tp{b-r^) = , 
J dHTppib) = 1, (4) 

and Ggg is the glue-glue interaction cross section. 

We are mostly interested in the initial transverse 
density distribution of the matter created at mid- 
rapidity in heavy-ion collisions. As the collision en- 
ergy increases, the production processes for this mat- 
ter involve gluons of smaller and smaller longitudi- 
nal momentum fraction x. As x decreases, the gluon 
density (gluon distribution function) grows [22]. At 
high gluon densities, gluons start to recombine [23], 
leading eventually to a saturation of the gluon den- 
sity [24]. For a Gaussian density distribution as in 
Eq. (1) gluon saturation becomes effective in the 
center of the proton first, reaching into the periph- 
ery only at still higher collision energies. This leads 
to a fiattening of the normalized density distribu- 
tion (1) near its center, accompanied by a growth of 
its rms radius. Gluon saturation effects can there- 
fore be mimicked by allowing the parameter B in 
Eq. (1) to grow with ^/s. This is what we suggest in 
the present article; it differs from the treatments in 
[14, 21] which held B constant. 

The inelastic nucleon-nucleon cross section tr™^' 
is obtained by integrating the inelastic gluon- 
scattering probability P{h) in Eq. (3) over all impact 
parameters: 



d%P{b) = d%(l- CXp[-aggTpp{b)] 



(5) 

In our model the y/s dependence of results from 
the combined energy dependences of agg and the 
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Gaussian width B of Tpp{b). This differs from the 
work [14] where the ^/s dependence of a'^^ arises 
entirely from the energy dependence of (jgg. 

With the exphcit form (4) for the binary gluon 
coUision density and the substitutions t = 6^/(4i?) 
and A = (Tgg/(47ri?), relation (5) can be written as 



AttB 



= /(A), 



where 



/(A) 



dt 



(l-cxp(-Ae-*)) 



(6) 



(7) 



is a monotonically increasing function with /(O) = 0. 

To completely specify our model we assume that 
A = (Tgg/(47ri?) is an energy independent constant, 
or (equivalently) that the Gaussian width B (which 
is proportional to the square of the proton radius) 
increases with collision energy proportionally to the 
inelastic nucleon-nucleon cross section Unn- This 
assumption makes intuitive sense: as the collision 
energy increases, wee partons in the tail of the pro- 
ton's density distribution get boosted to sufficient 
energy to contribute to the scattering between nu- 
cleons while a similar increase in the dense center of 
the proton gets tempered by gluon saturation. As a 
result, the proton gets bigger. The glue-glue inter- 
action cross section agg increases at the same rate; 
in fact, since P{b) is always smaller than 1 (due to 
the non- vanishing "punch-through" probability de- 
scribed by the last term in Eq. (3)), a-gg > cr™^' is 
always true.^ 

To fix A we use experimental input at fixed colli- 
sion energy that is sufficiently large for the Glauber 
model [2] and our gluon-saturation-based ideas for 
extrapolating to higher energies to work, but small 
enough that we can still make contact with the mea- 
sured and tabulated nuclear profiles at low ener- 
gies, without having to take into account energy- 
dependent swelling of the entire nucleus. 

At -ys = 23.5GcV, the inelastic nucleon-nucleon 
cross section has been measured as = 3.2 fm^ 



There is no fundamental reason for assuming A = 
agg / (AnB) to be constant. In fact, asymptotic freedom 
in QCD suggests that o-gg might grow more slowly with 
y/s than B, in which case A and /(A) would be decreasing 
functions of y's, and the proton size B would grow more 
strongly with energy than the inelastic nucleon-nucleon 
cross section. This would further increase the smoothing 
effects on initial state fluctuations at higher collision ener- 
gies to be discussed below. 



[25]. An extraction of the rms impact parame- 
ter t'rms = (6^)^^^ for inelastic proton-proton colli- 
sions from CERN ISR data has been reported by 
Amaldi and Schubert [26]: at ys = 23.5GeV they 
find forms = 1-03 fm. Our model allows to compute 
forms from the collision probability P{b) in Eq. (3) as 
follows: 



/rf^fofo^P(fo) 
Jd^bP{b) 

Jo dt(l-cxp(-Ae *)) 



(8) 



Equations (6) and (8) can be cast into the form 



B 



47r/(A) 45(A)' 



(9) 



which can be solved for both the energy-independent 
constant A and the proton width parameter B at 
V^=23.5GeV. We find A =1.62, corresponding to 
B(23.5 GeV) = (0.473 fm)^. At other collision ener- 
gies B can now be computed from 



(10) 



14.30 fm^ 

In Table I we collect a few representative values. 



TABLE I: The Gaussian width VB (Eq. (10)) for various 
collision energies. The values for ctnn a-t LHC energies 
(2.76 and 7TeV) were reported in [28, 29]. 



^^ (GeV) 


(mb) 


^/B (fm) 


23.5 


32 


0.473 


200 


42 


0.544 


2760 


62.1 


0.661 


7000 


71 


0.707 



Before describing how we distribute these "grow- 
ing nucleons" inside the colliding nuclei, let us sum- 
marize the novel features of our approach. In [14] 
ah of the energy dependence of crJJ^' on the l.h.s. 
of Eq. (5) was attributed to the factor agg in the 
exponent on the r.h.s. of Eq. (5). Here we split 
this energy dependence between the factors agg and 
Tppib) on the r.h.s. of Eq. (5), by letting both agg 
and the nuclcon size B grow with ^/s. Our specific 
model assumption (for which we have no convinc- 
ing argument other than simplicity) is that all three 
quantities (cJ^N' '^aa ^'^'^ ^) grow at exactly the 
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same rate. While this may not be completely cor- 
rect, some growth of the nucleon size B with y/s is 
suggested by theory [30] and required by experimen- 
tal data [26, 30, 31]: Amaldi and Schubert demon- 
strated a moderate growth of the proton inelastic 
interaction radius over the energy range probed by 
the ISR [26], and Fig. 19(E) in Ref. [30] together 
with the recent TOTEM result from elastic pp scat- 
tering at '/s = 7 TeV at the LHC provide experimen- 
tal evidence for a significant growth of the proton's 
elastic interaction radius by a factor ~ 2 between low 
ISR and top LHC energies, roughly comparable with 
the similar increase of a™^ over the same energy 
range documented in Table I. It is this growth of 
the proton radius that causes the smearing of event- 
by-event fluctuations in very high energy nuclear col- 
lisions to be discussed in the following. 



III. DISTRIBUTION OF NUCLEONS 

To describe the two incoming nuclei before the col- 
lision, we distribute the nucleons (1) with a distribu- 
tion pyi(ro) (where denotes the center of the nu- 
cleon) that is adjusted such that the resulting folded 
distribution [12] 



(11) 



reproduces the measured nuclear density profile for 
a nucleus of mass number A [27]. We take both pA 
and PA to be of Woods-Saxon form: 



PA{r) 



PAiro) 



Po 



Po 



exp 



ro-RA) 



(12) 
(13) 



The measured values for gold are [12, 27] 

po = 0.1695 fm"^ i?Au = 6.38fm, dAu = 0.535 fm. 

(14) 

We find that the choice po = po works with excellent 
accuracy (see Fig. 1). We perform a two-parameter 
fit for .Rau and d^n, using \/B = 0.473 fm in Eq. (1) 
{i.e. the value corresponding to the lowest collision 
energy in Table I, keeping in mind that the param- 
eters (14) are the result of low-energy scattering ex- 
periments). We thus find 



po = 0.1695 fm 



-3 



i?Au=6.42fm, dAu = 0.45fm. 

(15) 



These values are very close to those obtained in [12] 
for a box-like nucleon density profile. Fig. 1 shows 
the nucleon center distribution pAu (Eq. (13)) and 
the folded distribution pAu (Eq. (11)), and compares 
the latter to a Woods-Saxon profile (12) with the 
measured parameters (13). 
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FIG. 1: (Color online) Nuclear density pAu{r) (Eq. (11), 
solid line) and the distribution of nucleon centers /5Au(r) 
(Eqs. 12,14), dashed line), together with the Woods- 
Saxon distribution (12) with measured parameters (14) 
(dotted line). Please note that only the tail of the distri- 
bution is plotted, to better see the differences between 
PAu(r) and pAu(r). 



IV. ENERGY DEPENDENCE OF 
INITIAL-STATE FLUCTUATIONS IN 
AU+AU COLLISIONS 

The smooth dashed curve in Fig. 1 can be inter- 
preted as the (unnormalized) probability (i.e. the 
squared modulus of the single- nucleon wave func- 
tion) for finding a nucleon centered at position r in- 
side a gold nucleus. At high energies, where the time 
for the colliding nuclei to pass through each other is 
very short, each nucleus-nucleus collision can be see 
as performing an instantaneous position measure- 
ment of the positions of all struck nucleons inside 
the nuclei (i.e. of the A-particle wave function in 
position space, integrated over the positions of those 
nucleons that pass through unscathed). With few 
exceptions (see e.g. [32]), Monte Carlo implemen- 
tations of these ideas assume uncorrelated nucleons 
(except for, in some cases, an excluded volume corre- 
sponding to an infinitely strong hard-core repulsion 
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at short distances), i.e. a factorization of the A- 
nucleon probabihty density into single-nucleon prob- 
abihtics, given by the normahzcd version of Eq. (13). 
This allows to calculate the probabilities for dis- 
tributing the A nucleons inside each nucleus inde- 
pendently. 
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FIG. 2: (Color online) The transverse entropy den- 
sity s{r±,To =0.4fm/c) (in fm~'^) at mid-rapidity on an 
8fm X 8fm grid from the MC-KLN model, for a cen- 
tral {b = 0) Au+ An collision. The top graph uses disk- 
like nucleons (Eq. (20) with a'^^ = 42mh), whereas the 
middle and bottom plots use Gaussian nucleon thickness 
functions with y/B = 0.54A and 0.707 fm, corresponding 
to collision energies of — 200 GeV and 7 TeV, respec- 
tively (see Table I). All three plots use the same collision 
event, i.e. identical spatial distributions for the nucleons 
in the two nuclei. 

We here use the MC-KLN model described in 



Ref. [11]. It calculates initial gluon production by 
folding the unintegrated gluon distributions of the 
two nuclei with a cross section for gluon fusion [33] 
(see Refs. [4, 11, 12, 15] for details). The unin- 
tegrated gluon distributions depend on transverse 
position through the gluon saturation momentum 
[11, 12] 



0.01 



where 



TAiri_) = J dzpA{r^,z) 



(16) 



(17) 



is the nuclear thickness function of the correspond- 
ing nucleus at the same position, Ta^q = l-53fm~^, 
and X = pt exp(±y) / y/s are the longitudinal momen- 
tum fractions of the gluons from the two colliding 
nuclei that fuse into a produced gluon with rapid- 
ity y and transverse momentum px. The Monte 
Carlo sampling of the position distribution pAifo) in 
Eq. (11) leads to event-by-event fluctuations in the 
nuclear thickness function TA(r±) that (via Q'^{r±)) 
are reflected in event- by- event fluctuations of the 
produced transverse gluon density which, after ther- 
malization, generates the initial entropy and en- 
ergy density profiles for the hydrodynamic evolution. 
Through Eq. (11), the shape and character of these 
fluctuations is affected by the nucleon density profile 
Ppir)- 

Integrating Eq. (11) over the longitudinal position 
z yields 

TAir±)^ [ d\±ofAir±o)Tp{r^-r^o), (18) 



which, after distributing the nucleons in the nuclei 
by Monte Carlo sampling the distribution ^^(ro), 
becomes 



TA{r^)=J2Tpi\r^-n^\)- 



(19) 



Here the terms under the sum are Gaussians (2) with 
a width B that increases with as described in 
Sec. II. Previous implementations of the MC-KLN 
model [11; 12] instead used a disk-like nucleon thick- 
ness function; 



^incl 
'NN 



-inol 



(20) 



In Figs. 2 and 3 we compare the initial transverse 
gluon density distributions for a central Au-|-Au col- 
lision with disk-like and Gaussian proton thickness 



5 




-8-6-4-2 2 4 6 8 -8-6-4-2 2 4 6 8 
X (fm) X (fm) 

FIG. 3: (Color online) Influence of the growth of the proton with collision energy on the initial transverse entropy 
density s(rx, ro = 0.4 fm/c) (in fm^'') for a typical central (b — O) Au+Au collision from the MC-KLN model. All 
four panels use identical nucleon distributions in the two colliding nuclei, but the width ^/B of the nucleon thickness 
function changes with ^/s (indicated above each panel) as listed in Table I. 



functions and for different collision energies, taking 
into account the growth of the width ^/B of the nu- 
cleon thickness function with increasing cohision en- 
ergy. The upper two panels in Fig. 2 illustrate the 
significant smoothing of the initial fireball density 
distribution at small length scales that is caused by 
replacing the disk-like nucleons with Gaussians. The 
small but pervasive discontinuities in the density 
profile using disk-like nucleons are unphysical but 



constitute a major technical stumbling block for vis- 
cous hydrodynamics which cannot run stably with 
such discontinuous initial profiles (unless an addi- 
tional smoothing step, involving another unphysical 
parameter, is applied first). With Gaussian nucleons 
these discontinuities arc naturally washed out. 

However, the initial density distributions shown 
in the middle panel of Fig. 2 and in the upper two 
panels of Fig. 3 still feature large (physical) inhomo- 
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geneities and density gradients which, for very early 
starting times of the hydrodynamic evolution (where 
their effect on the viscous pressure is enhanced by 
a factor l/r [34] that arises from the longitudinal 
Bjorken expansion), can drive viscous hydrodynam- 
ics outside its region of applicability. The length 
scale of these initial state density fluctuations thus 
generates a lower limit for the hydrodynamic start- 
ing time tq. At earlier times, even if the system 
started out in a state of local equilibrium, viscous 
forces would drive the system so far away from ther- 
mal equilibrium that the macroscopic hydrodynamic 
framework breaks down. 

The lower two panels in Fig. 2, as well as the pro- 
files shown in Fig. 3, demonstrate that the grow- 
ing nucleon size smoothes the initial state density 
fluctuations and significantly increases the length 
scale over which the initial pressure profile fiuctu- 
ates. This implies that for heavy-ion collisions at the 
LHC viscous hydrodynamics can be applied starting 
from significantly earlier times than at RHIC ener- 
gies. 

V. FURTHER DISCUSSION AND 
CONCLUSIONS 




FIG. 4; (Color online) The ratio 

(e„(200GeV))/{e„(2.76TeV)) of the ensemble-averaged 
harmonic eccentricity coefficients (e„) (n = 2, . . . , 15) 
at y/s = 0.2 and 2.76 A TeV, for Au-|-Au collisions with 
MC-KLN initial conditions as a function of the number 
of participating (struck) nucleons, A'part. The plot is 
based on 95,000 events, binned into 10 equal size A^part 
bins containing 9,500 events each. Smoother curves 
require higher statistics. 



The smoothing effects of a larger effective nucleon 
size at higher energies also influence the anisotropic 
flow generated in heavy ion collisions. Initial-state 
density fluctuations entail event- by- event shape fluc- 
tuations for the initially produced fireball. As dis- 
cussed and studied in detail in Refs. [19, 35-37], 
these can be characterized by a set of harmonic ec- 
centricity coefficients e„ (n = l,2,...) which drive 
higher order harmonic fiow coefficients w„. 

Figure 4 shows the ratio of the ensemble averages 
(£„) of some of the first 15 harmonic eccentricity 
coefficients at two different energies. The e„ are de- 
fined in terms of the initial transverse energy density 
profile e(rx,0) through [35, 38] 



J.2 gm0 



e(rj_. 



J cPr± r\ e(rj^. 



(21) 



Fig. 4 is for Au-|-Au colhsions at RHIC and LHC 
energies as a function of collision centrality, using 
the MC-KLN model with Gaussian nucleons whose 
width grows with energy. We see that, except for the 
ellipticity £2, all higher order harmonic eccentricity 
coefficients are larger at RHIC than at LHC, with 
the ratio between the two energies increasing with 
harmonic order n to about 1.5—1.6 for 71^=15. This 
is a straightforward reflection of the higher degree 



of granularity in the initial density profiles at RHIC 
than at LHC, caused by the smaller effective nucleon 
size. The ellipticity £2 bucks this trend: in contrast 
to the higher order harmonics which are driven by 
fluctuations, the ellipticity is dominated in all but 
the most central collisions by the elliptic geometric 
deformation of the nuclear overlap zone. The geo- 
metric contribution to the ellipticity actually grows 
with collision energy, due to the slight swelling of 
the entire nucleus caused by the growth of its con- 
stituent nucleons. Consequences for the extraction 
of the QGP shear viscosity from anisotropic flow co- 
efficients will be explored in future work. 

In summary, when taking into account the 
swelling of the nucleons at higher collision energy 
caused by gluon saturation effects, we expect the el- 
liptic fiow to slightly increase, but the higher order 
flow harmonics to decrease at the LHC compared to 
previous predictions. This effect increases with the 
harmonic order of the flow. Furthermore, smoother 
event-by-event initial conditions extend the range of 
validity of a viscous hydrodynamic description to 
earlier times, so taking into account the swelling of 
the nucleon from RHIC to LHC energies further adds 
to the applicability of viscous hydrodynamics for 
the description of the dynamical evolution of ultra- 
relativistic heavy-ion collisions. 
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